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Abstract. In this study, with the aid of N-Body simulations based on quasi-cosmological initial conditions, we have 
followed the formation and evolution of two early-type galaxies of different total mass, from their separation from 
global expansion of the Universe to their collapse to virialized structures, the formation of stars and subsequent 
nearly passive evolution. Using the PD-TSPH N-Body code, we developed the simulations of two model galaxies 
in quasi-cosmological context. The cosmological background we have considered is the Standard CDM. The 
models are made of Dark and Baryonic Matter in standard cosmological proportions (9:1), and have significantly 
different initial total mass, i.e. 1.62 x 1O 12 M (A) and 0.03 x 1O 12 M (B). Particular care has been paid to the 
' star formation process, heating and cooling of gas, and chemical enrichment. Star formation is completed within 

, the first 3 Gyr in Model A, whereas it lasts longer up to about 4 Gyr in Model B. The models are followed for a 

long period of time, i.e. 13 Gyr (Model A) and 5 Gyr (Model B): in any case, well beyond the stages of active star 
formation. The structural properties of the present-day models are in good agreement with current observations. 
The chemical properties, mean metallicity and metallicity gradients also agree with available observational data. 
Finally, conspicuous galactic winds are found to occur. The models conform to the so-called revised monolithic 
scheme, because mergers of substructures have occurred very early in the galaxy life. Our results agree with those 
obtained in other similar recent studies, thus strengthening the idea that the revised monolithic scheme is the 
|' right trail to follow in the forest of galaxy formation and evolution. 
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.' r , , „ , , . . Il975l lArimoto k Yoshiil 119871 Bressan et alJll994l) ever 

K> (Jne oi the maior challenges m modern astrophysics is ,, „ , , . , r r < i • i 

/\ ... . . since followed by quiescence. In favour ol this scheme 

S_i ' to understand the origin and evolution oi galaxies, the ., , ; . , , . . , , , 

3 . ..... . . T , . are the observational data that convincingly hint lor old 

bright elliptical ones m particular. In a universe domi- , , ... , , . , tptT- — 1| nnnr J r 

° , . z, , , _ . , , f^n^-K ,\ i i an d homogeneous stellar populations (see Chiosi 2000, for 

nated by Cold Dark Matter CDM) and, as recent obser- . „ , , , . , N T , . . 

a recent revie w ol the su bject), ft is worth mentioning, 

vations seem to suggest, by some kind of Dark Energy , , , , rr — "Z — "^1 l. nno i. , — , 

. . „ A i • i a i however, that Kaunmann et al. ( 1993) and Barger et al. 

m the form oi a non-zero Cosmological Constant A, and 177^71. r , i r i , • • , i 

. , ° (1999) argue for the presence ol recent evolution m the 

containing a suitable mix oi baryons and photons, cos- . „ , ,. n ... , , . 

° _ iii • • i n stellar populations oi elliptical galaxies, 

mic structures are formed by the gravitational collapse _., , . , . . . . . . 

, .,. ,• i r 21 he hierarchical scenario suggests instead that 

ol Dark Matter and are organized m a hierarchy ol com- v / , . , °° . 

. /i i \ • • i i-ii i- • . massive galaxies arc the end product ol subsequent merg- 

plexes (halos) msidc which baryons dissipate their energy . , . , , , , 

... . c i • T ii • a. 1 ers ol smaller subumts, over time scales almost equal to the 

and collapse to form luminous systems. In this context, ..... , . 



modelling the formation of early-type galaxies with simu- 
lations taking into account the dynamics of Dark Matter 



Hubble time in the particular cosmological model used; as 
the look back time increases, the density in comoving space 

. . ... ,. ,. . of bright (massive) elliptical galaxies should decrease by 

halos and gas, radiative cooling ol baryons, star formation, r , ° , „ , |„ r , '., — „ " — ik j rr. — ~ — ■ ; — ri 

. , , , . , , , , . P „ . a factor 2 to 3 e.g. White & Rccs 1978: Kaunmann et al. 

and gas loss by galactic winds can be reduced to following i. - _ r r ~Tr ! — — , — ~. ; . 

schemes JSQ Eo02 : Schade eTallB . ^ ^° Ur ° f thlS ™ " e SOme obsei ;yational evi- 

/i\ m i-,7 ■ • r i r dcnccs that the merger rate likely increases with ~ ( 1 +z) 

(1 the monolithic scenario oi galaxy formation sup- I,, , , — ■ — ni. ~jt , .. .., , . . „ , 

. ,. , ,, , ,i , , , • r , (Patton et al. 1997 together with some hint for a color- 
poses and predicts that all early-type galaxies form at . . , . ° _ . „_ , . ,. . , 
... , . , . , n , i . , . structure relationship for E & SO galaxies: the color be- 
high rcdshiit by rapid collapse and undergo a single promi- . . , . . . ° „ . 

or- comes bluer at increasing complexity ol a galaxy struc- 

Send offprint requests to: E. Merlin ture. This could indicate some star formation associated 



2 



E. Merlin & C. Chiosi: Formation and evolution of early-type galaxies 



to the merger event. Finally, there are the many successful 
numerical simul ations of galaxy encounter s, mergers and 
interactions fe.g. lBarnes fc Hernauistll996t) . Nevertheless, 
contrary to the expectation from this model, the number 
density of elliptical does not seem to decre ase with the 
redshift, at least up to z ~ 1 l|lm et al.lll996|) . 

(3) There is a third scheme named dry merger view, in 
which bright ellipticals form by encounters of qu iescent, no 
star fo rming galaxies. This view is advocated bv lBell et alJ 
( 2004^ who finds that the -B-band luminosity density of the 
red peak in the colour distribution of galaxies shows mild 
evolution starting from z ~ 1. As old stellar populations 
would fade by a factor 2 or 3 in this time interval, and the 
red colour of the peak tells us that new stars are not being 
formed in old galaxies, he argues that this mild evolution 
hints for a growth in the stellar mass of the red sequence, 
either coming from the blue-peak galaxies in which star 
formation is truncated by some physical process, or by 
"dry mergers" of smaller red, gas-poor galaxies. 

(4) Finally, a fourth hybrid sc enario named revise d 
monolithic has been proposed by ISchade et al.l l|l999l) . 
who suggests that a great deal of the stars in massive 
galaxies are formed very early-on at high z and the remain- 
ing few ones at lower z. The revised monolithic ought to 
be preferred to the classical monolithic, as some evidences 
of star formation at 0.2 < z < 2 can be inferred from the 
presence of the emission line of [Oil] , and also as the num- 
ber frequency of early-type galaxies up to z ~ 1 seems to 
be nearly constant. 

However, a sharp distinction among the various sce- 
narios could not exi st in reality. As pointed out by 
lLonghetti et al.l l|200(l) early- type galaxies in isolation and 
in interaction (such as pair and shell galaxies) share the 
same distribution in diagnostic planes such as HP vs 
[MgFe], the classical tool to infer age differences. This 
means that secondary episodes of star formation may not 
only occur in merging (dynamically interacting) galaxies, 
but also in the isolated ones because of internal processes, 
thus supporting the Schade revised monolithic scheme. 
Captures (mergers) of small satellites by a galaxy born 
in isolation according to the monolithic scheme are rea- 
sonably possible (e.g. the Milky Way which is currently 
capturing the Sagittariu s dwar f galax y ). In any case, as 



of the constituent stellar populations. It can be shown 
that the colour-magnitude relation stems from a mass- 
metall icitv sequence and no t an a ge se quence, see for in- 
stance iKodama fc Arimotol <)l997l) and IChiosi fc Carrarol 



Chiosi fe Carraro 



Tantalo fc Chios 



thorou g hly discussed bvlchiosil l200o|) , 
<|2002h iTantalo fc Chiosil (l2004af and 
Ij2004bl) . the age and intensity of the last episode of star 
formation (measured by the fractionary mass engaged in 
newly formed stars) cannot exceed some stringent limits, 
a few per cent and about one third of the Hubble time 
even in a massive elliptical, otherwise the typical broad 
band colours cannot be matched; they would be too blue 
compared to observations. 

In addition to all this, elliptical galaxie s are known 
to ob ey to a colour-magnitude relation l|Bower et alJ 
119921) : the colours become redder at increasing luminos- 
ity (and hence mass) of galaxies. The relation is tight 
for cluster galaxies and more disperse for field galax- 
ies. The colour of a galaxy simply mirrors the colours 



( 2002) and re ferences therein . 

Long ago, iLarsonl l|l975|) advocated the supernova 
driven galactic wind model to explain the colour- 
magnitude relation. He argued that massive galaxies re- 
tain gas and form stars for longer periods of time, and 
become richer in metals than the low-mass ones. In any 
case the star forming period even in the most massive 
galaxies seemed to be confined at very early epochs (say 
within the first Gyr to be generous). This kind of model 
has been extensively used to describe and predict the 
chemo-s pectro-photometric p r operties of elliptical galax- 
ies ('e.g. iRTelsan et alJ ITflflfl [S ITantalo et alJll 
1993 IChiosi et alJll998l) 



The short duration of the dom- 
inant star t o rmmg period has also been suggested by 
iBower et alJ l)l992i) on the basis of very simple, empiri- 
cal arguments: the vast majority of stars are formed at 
z > 2 and subsequent activity if any should have occurred 
at z > 0.5. 

Over the ye ars, the situat ions became more intrigued 
because (i) the ILarsonl ljl975|) model got in contrast with 
the observational constraint that massive elliptical have 
on the average [a/Fe] ratios higher (super-solar) than the 
low mass ones (solar or below solar), thus implying the 
opposite trend for the duration of the star forming period 
as a mply discussed bylchiosi et all l|l998|) . Ichiosil (|2000^ 
and IChiosi & Carrarol l|2002|k fii) of the results from ab- 
sorption line indices, which seem to indic ate that large 
age d iff erence could e xist. According to iBressan et alJ 
(Il996h. ITantalo et all ill 998lL ITantalo fc Chiosil (l2002[) . 
ITantalo fc Chiosil l|2004a|) and ITantalo fc Chio'sTl|2004l| . 
the following picture emerges: all elliptical galaxies have 
been formed far back in the past, but have undergone dif- 
ferent histories of star formation. Massive ellipticals had 
a single burst-like episode of stellar activity ever since fol- 
lowed by passive evolution, whereas the low mass ones 
stretched their star formation history over long periods of 
time in a series of bursts. This view has been supported 
by the N-Body simulations of IChiosi fc Carrarol ( 20021) 
who find that the duration and mode of star formation 
is driven by the total mass (Dark plus Baryonic Matter) 
and/or the initial density of th e galaxy. In addition to 
this the IChiosi fc Carrarol l|2002f > models while predicting 
copious galactic w inds much im proved upon the simple- 
minded scheme hv ILarsonl ill 97,^ and ruled out the point 
of inconsistency with the abundance ratios. We will came 
back to these issue later on. It is also supported by the 
most recent Spitzer-data (and SIRTF) in the far infrared, 
which have revealed the existence of very massive galax- 
ies already in place at redshift z ~ 6; in addition to data 
HST observations hav e also brought into ev idence galaxies 
in place at z ~ 7 - 8 llBouwens et al.ll20o4 - 

All this can be understood within a model ruled by 
violent relaxation taking place in the remote past thanks 
to which clusters of just born stars can merge together 
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into even denser clusters fe.g. IMacchetto et al.lll99(j) . In 
this context it is als o worth mentioning that according to 
IBundv et all l|2005(l (i) the mass of star forming galaxies 
increases with redshift ( "downsizing") thus strongly weak- 
ening the hierarchical scenario as the main mechanism for 
forming giant ellipticals; and (ii) even if "dry mergers" 
clearly occur l|van they cannot play an 

important role in the assembling history of massive quies- 
cent galaxies as the expected number density of massive 
galaxies at low redshift does not agree with observational 
data. Furthermore the mass-downsizing of star forming 
galaxies seems to be almost independent of the environ- 
mental density, whereas one would expect a strong density 
dependence in the dry merger rate. Therefore, either clas- 
sical or revised monolithic scenarios seem t o be the right 
frame for the galax y formation mechanism llPeebleJ 2002: 
ISchade et allll999l) . 

B asing on the pione er studies by iKat d l)l99l[) 
and iKatz fc Gunnl l|l99ll) . a nd later by many au - 
thors among who m we r ecall IKawatal Jl999l l2001albh : 
iKawata & Gibson! l|2003h and iKobavashil ^POSlT lve 
present nume rical simulations ca rried out with the PD- 
TSPH code ijCarraro et alJll99^) in order to investigate 
the formation mechanism on galactic scales. The scenario 
we have in mind and we intend to prove with the present 
simulations is the following one. At a certain (high) red- 
shift a massive perturbation made of Dark Matter and 
baryons detaches from the Hubble flow and collapses on 
its own. It becomes a massive proto-galaxy, rich of sub- 
structures which moving inside the common gravitational 
potential well merge, form stars and eventually give rise 
to a single entity (the galaxy). The process is complete 
at high redshift (say before z ~ 2) and from now on 
the galaxy evolves in "isolation"; subsequent captures of 
small satellites are possible without significantly altering 
the overall structure and evolution. In the case of a low- 
mass perturbation the stellar activity is prolonged over 
long time scales. Mergers between two galaxies may occur 
but they are a sort of rare, spectacular event and not the 
basic mechanism for assembling massive ellipticals. 

The plan of the paper is as follows. Section [5] presents 
the procedures we have used to derive the initial condi- 
tions for our models; the procedure we have adopted is 
outlined and compared with other techniques available in 
literature. Section |3| describes some assumptions we have 
made for some fundamental physical processes that are in- 
cluded in the numerical simulations, namely cooling, star 
formation and initial mass functions, heating and chem- 
ical enrichment, and shortly outlines the numerical im- 
plementation of these processes in the N-Body PD-TSPH 
code developed by the Padova group over the years, to- 
gether with a brief description of the probabilistic view of 
star formati on, heating, cooli ng and chemical enrichment 
proposed bv lLia et alJ (|200^ . Section 0] presents the key 
parameters of the models we plan to calculate, and in par- 
ticular it deals with the initial conditions derived from the 
method described in Section [5] Section presents the re- 
sults: the spatial evolution, the star formation history, the 



energy balance, the chemical enrichment and other details. 
Finally, Section contains some concluding remarks. 



2. Search of the initial conditions 

The choice of the initial conditions is perhaps the fun- 
damental step on which much of the future evolution of 
struct ures (galaxies) will depend. As noted bv lPower et alJ 
( 2003), despite the popularity of cosmological N-body sim- 
ulations in the last years, there is very little detail in lit- 
erature about setting up initial conditions. 

The scheme we have adopted can be defined as "quasi- 
cosmological" . It starts from realistic simulations of the 
cosmological evolution of the Universe that are suitably 
adapted to the numerical simulations of galaxy formation 
and evolution. 



2.1. Quasi-cosmological initial conditions 

Our i ni tial conditions stem from the stud i es of iKatzl 
l|l99ll) . IKatz Gunrl lll99lh and IKawatal l)l999h . We 
adopt the Standard Cold Dark Matter Universe with 
Q(z) = r2o, the matter density at the present time. 
Therefore, the basic parameters of the model are Hq, 
Oo, and eg (the rms mass fluctuations of the present-day 
Universe that normalizes the density fluctuations and fix 
the initial redshift of the simulations via the amplitude of 
the density fluctuations) 1 . These models do not consider 
the presence of the Cosmolo gical Constant A as i ndicated 
by the recent W-MAP data ((Bennett et aljEoO^ . 

2.1.1. Preparing the quasi-cosmological models 

To derive the initial conditions (positions and velocity 
field) of Dark and Baryonic Matter particles of the pro- 
galaxy, we start from simulating an initial large grid rep- 
re senting the Universe using the public software Grafic2 
bv lBertschingerl l|200lh . The initial Universe grid is a very 
large cube with ~ 40 comoving Mpc per side. In the grid 
we select a peak of over-density which is assumed as the 
centre (or very close to it) of a new cubic sub-grid of 
smaller size. In practice, we fix the size of the sub-grid 
looking for the distance from the peak at which the over- 
density goes to zero. This distance is taken as half the size 
of the cube dimension I. In this sub-grid both Dark and 
Baryonic Matter are present. Their respective barycentres 
do not necessarily coincide each other nor are both coinci- 
dent with the geometrical centre of the cube. Then we iso- 
late this cubic volume, and using Grafic2 multi-level ini- 
tial conditions generator we improve upon its resolution, 
thus describing small-scale perturbations in great detail, 
superposed to the large scale perturbations described in 
the first-level grid. Subsequently, using the radial distance 
of each particle (Dark Matter and Baryons) from the cen- 
tre of the cube, we select the sample of particles contained 



1 Al l equations used in this section are from iPadmanabhanl 
(1993) unless otherwise specified. 
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in a sphere with radius equal to 1/2. The volume of the 
sphere is about half of the volume of the cube and ac- 
cordingly the number of particles in the sphere is about 
half of the number of particles in the cube. It is worth re- 
calling that, due to the cosmological perturbations, as the 
baryons in the cube could have a slight offset with respect 
to Dark Matter, the same could happen with the spheres 
of baryons and Dark matter. Therefore, the numbers of 
particles of the two matter species can be slightly differ- 
ent, and not be exactly concentric. This however is less of a 
problem because first it is physically understandable (per- 
turbations of the density field and infall of baryons in the 
Dark Matter potential wells), and second the difference in 
the numbers of particles is always small. By replacing the 
cube with a sphere we avoid spurious tidal effects on the 
particles. The sphere represents our proto-galaxy in very 
early stages. It is worth noticing that centering the new 
coordinate system at the geometrical centre of the cube 
rather than at the mass-weighed centre leads to the pos- 
sibility that the peak of the density perturbation does not 
coincide with the centre of the sphere. The particle posi- 
tions are then referred to the new system: the distances 
provided by Grafic2 in comoving Mpc are translated to 
proper distances in Mpc according to 



d>pr Qstart ^ d c 



(1) 



where a sta rt = (1 + Zstart)^ 1 is the value of the expansion 
factor at the initial redshift, d comov is the distance in the 
comoving frame of reference and d pr is the proper distance. 
The velocities are then referred to the centre of the sphere, 
dropping the motion of the sphere as a whole and leaving 
to each particle only the peculiar motion. To each particle 
we add the velocity of the Hubble flow, which is simply 
given by 



Vflow — H(z s tart) x d 



pr 



(2) 



(this velocity is radially oriented). Note that in order to 
calculate H(z) one has to assume a cosmological model. 
In the Standard CDM cosmology 

H(z) = H ^/n m (l + z) 3 . (3) 

Finally we add a certain amount of rigid rotation. The 
associated spin parameter A is 



A = 



J\E\ 1 ' 2 



(4) 



GM 5 / 2 

where J is the angular momentum, E is the initial bind- 
ing energy, and M the total mass of the system. Typical 
values of A range from 0.02 and 0.08 (|Whitdll98^ . which 
corresponds to angular velocities of the order of fractions 
of a complete rotation over time-scales as l ong as about 
ten free-fall time-scales l)Carraro et alJ ll998). As the free- 
fall time-scale depends on the density according to 



1 



it strongly depends also on the initial redshift. The back- 
ground density for matter scales indeed as (1 + z) 3 . 
Introducing the rigid rotation somehow simulates and 
takes into account the effects of fluctuations with wave- 
length longer than the typical size of the proto-galaxy (e.g. 
lKawatalll999|) : neglecting rotation would imply that the 
matter surrounding the proto-galaxy has no effect on its 
evolution. 

2.1.2. Densities and masses: some cautionary remarks 

The theory of spherical collapse (cfr. IPadmanabharll 993) 
would allow us to estimate the approximate redshift at 
which turn-around and collapse of an homogeneous and 
expanding sphere should occur under the action of self- 
gravitation. The relations of interest here are 



1 + Z turn = 0.57(1 + Zi)Si 

1 + z coU = 0.36(1 + Zi)8i 



(6) 



(5) 



where Si is the mean over-density of the sphere (proto- 
galaxy) at the initial redshift. But these values are mere 
approximations of reality. 

First of all, the classical theory of spherical collapse 
stands on top-hat perturbations, in which the particle ve- 
locity field is determined only by the Hubble flow. The 
presence of peculiar velocities and density perturbations 
strongly alters the predictions of this simple theory. If the 
initial conditions are such that there is a single, strong 
peak in the density field, a more reliable model could 
be th e secondary infall model (e.g. iMiller fc Bertschingerl 
1986J). Models for the spherical collapse taking the initial 
kinetic energy of peculiar mot i ons in to account do actually 
exist, see e.g. IPadmanabhanl l)l993h . As these models go 
beyond our aims, they will not be consider here. Second, 
cutting the sphere out of the cosmological Universe-grid 
and isolating it in a void space implies that the mean den- 
sity of the Universe coincides with the mean density of the 
sphere. Therefore, assigning a certain over-density to the 
sphere will force it to follow an evolutionary path typical 
of a Universe with £1 >1. This means that it will col- 
lapse on a faster time-scale. In cosmological simulations, 
the problem is nowadays solved using periodic boundary 
conditions. To cope with this point of uncertainty in our 
models, we have limited ourselves to consider only low 
initial over-densities (< 0.1). 

Furthermore, even if an univocal relationship should 
exist between redshift and proto-galaxy over-density, in 
reality fluctuations around the mean value of 6 are pos- 
sible, the fluctuations being larger at decreasing mass of 
the proto-galaxy. 

Finally, the spherical symmetry and sharp confinement 
of matter in the sphere we have adopted to avoid complica- 
tions due to possible irregular shapes of the proto-galaxy, 
imply that particles in the outskirts of the galaxy will meet 
sudden vanishing of density and velocity at the surface of 
the sphere. 
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Given all these uncertainties, it is clear that the results 
of the simulations will be only loosely related to theoretical 
predictions such as expected collapse redshifts. 

Considering the total mass as a free parameter would 
allow to investigate a wider range of possibilities, but 
such a choice would somehow weaken the above quasi- 
cosmological scheme and the peculiar velocities (which 
depend on the overdensity) should be properly adjusted. 
Therefore we prefer to assign the sphere the mass given 
by the "correct" value of the mean overdensity output by 
Grafic2. 

2.2. Comparison with other initial conditions 

It is worth commenting here the main diffe rence between 
our in itial c onditions and thos e adop ted bv lCarraro et alJ 
l|l998|) and IChiosi fc Carrarol l)2002|) . They start with a 
spherical model containing Dark and Baryonic Matter 
with masses in the standard proportions 9:1. They set the 
initial density profile for the Dark Matter component to 
be 



P( r ) = Pc — 



(7) 



where r c and p c are radius and den sity of a small centra l 
sphere. It resembles the well-known iNavarro et all l)l996|) 
law (see below) in the central regions. This indeed, in- 
dependently of cosmological models, initial spectrum of 
perturbations and total mass, this profile is oc r _1 in 
the central regions and oc r~ 3 in the external ones. The 
spatial positions of Dark Matter particles are obtained 
by MonteCarlo deviations from the distribution law. The 
initial velocities of Dark Matter particles are obtained 
from the velocity dispersion cr(r) assuming equi-partition 
among the three components, i.e. 



v(r) = \a{r) = Jp c r c Gr ln(^) (8) 



which is obtained by inserting eqn. Q in 



p{r)a{rf 



R - GM(r') , 



(9) 



l)Binnev fc Tremaindll987]) . 

The gas particles (baryons) are then homogeneously 
and randomly distributed in the spherical halo of cold 
Dark Matter with null velocity field to simulate the infall 
of primordial gas int o the potential well of Dark Matter 
l|White fc Reeslll978|) . 

In practice, this is equivalent to start with a halo of 
Dark Matter already detached from the Hubble flow which 
begins to collapse carrying along the baryonic matter. 
Even if all this sounds physically reasonable, these initial 
conditions already contain the solution of the problem: 
a self-gravitating, collapsing halo of Dark Matter whose 
density profile already resembles the final structure one is 
looking at. Moreover, in the CDM cosmology catching up 
of baryons by the potential wells of Dark Matter happens 



much earlier than the time at which a proto-galactic halo 
detaches from the Hubble flow. 

The initial radius is set to be equal to the virial radius 
i?2oo of the proto-galaxies, which is defined as the radius 
of a sphere with mean density equal to 200 times the mean 
background density of the Universe, p u , which in the case 
of a Standard CDM model (f^o = 1) is given by 



Pu(z) 



ml 

8nG 



(1 + z) 3 



(10) 



with obvious meaning of the symbols, since the min- 
imum density requi r ed for collapse is a bout 200p u (z) 
i Navarro et aDll996t |Padmanabhan]ll993j) . Finally, star 
formation is supposed to occur only after the virialization 
of the CDM halo. 

To summarize, the present initial conditions differ from 
those used in past similar studies carried out by the Padua 
group with the same software (e.g. lChiosi fc Carrarol2 002) 
at least in the three aspects: 

— The initial velocities taking into account the peculiar 
motions and the Hubble flow provide a more realistic 
description of the complexity of the galaxy formation 
process as compared to the "static" model. 

— The variation of the perturbation amplitude and spec- 
trum allows us to better explore their effects on the 
resulting galaxy models and properties. 

— The peculiar velocities are set to their initial value in a 
self-consistent way securing the correct g rowth of th e 
perturbations on all scales of interest (e.g. lKat3ll99lj) . 

Modern numerical cosmological simulations include 
more sophisticated methods, such as evolution in comov- 
ing space and resimulation at different degrees of accu- 
racy and resolution (see, e.g. lKawata fc Gibsonl200^) . The 
present scheme has to be intended as the initial step to- 
wards future improvements. Nevertheless we are confident 
that the present method, despite its limitations, is reason- 
ably suited to study the formation and the morphological, 
dynamical and chemical evolution of early-type galaxies as 
isolated objects. 

3. Basic physics in galaxy models 

In this section we shortly summarize the basic phys- 
ical processes that are thought to lead galactic for- 
mation and evolution, together with some details on 
their implementation in the fully lagrangia n PD-TSPH. 



Carraro et alJ 



2002) and 



More de tails on the cod e can b e found in 
lll998l). iBuonomo et all l|2000|) . iLia et al 
IChiosi fc Carrarol l)2002fi . 

Key physical ingredients ingredients of the whole prob- 
lem are the gravitational force that acts on all kinds of 
particle and is imple mented by the classical Tree-code of 
iBarnes fc Hutl l)l986j) with quadrupole expansion, and the 
hydrodynamic processes, which instead act only on the 
energy and motions of the gas particles together with star 
formation, chemical evolution, gas cooling by radiative 
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processes, and finally gas heating by energy injection from 
supernova explosions, stellar wind, UV cosmic radiation 
and others. 

The system of gas particles can be viewed as a colli- 
sional fluid governed by the classical equations, which in 
lagrangian notation are 



processes than simple shock waves. There are many pos- 
sible formulations fo r this additional term, among which 
we adopt the one by iMonaehan fc Lattanziol l)l985|) 



n 



Ocpij Cij 



where 



Pij 



Jij 



(18) 



dp 

~db 
dv 

a 

du 
~dl 
V 2 $ 



— pV v 

--VP - v$ 

p 

Vu 

P 

4nGp 



S 



(11) 
(12) 

(13) 
(14) 



i.e. the continuity or mass conservation equation, the 
Euler or momentum conservation equation, the energy 
conservation equation, and the Poisson equation linking 
the gravitational potential to the medium density, respec- 
tively. The companion equation of state is that of a perfect 
gasP=^. 

Gas-dynamics is described w ith t he SPH formalism 
bv IMonaehan fc Lattanziol l|l985f ) andlHernauist fc Katzl 
ill 9891). as im ple mente d bv ICarraro et alJ (jl998^ 7 
Buonomo et~aT] jExxj), iLia et alJ l)2002|) and 

Chiosi fc Carrara l|2002|) : in particular, a probabilis- 
tic approach to non-adiabatic processes such as star 
formation and gas restitution is used (see below). The 
term S in ean. H13|) is the energy source (otherwise known 
as source function) from all heating and cooling processes. 
It is written as 



S = 



A-r 
p 



(15) 



where A and T are the cooling and heating rates, respec- 
tively (see 13. II and 13. 31 for details). 

To derive the mean value of any physical quan- 
tity f(r) over a certain interval we adopt the so- 
called "gather/sca t ter" functional relation suggested by 
iHernauist fc Katzl l)l989j) : 



<f(r) >= J f(r')±{W(r-r\h(r'))+W(r-r',h(r))}dr', 

(16) 

where h is the smoothing length, which gives the exten- 
sion of the volumes over which a suitable space average is 
mad e. For the kernel W(r) we a dopt the form proposed 
bv lMonaghan fc Lattanziol <|l985|) 



W(r, h) 



1 

irh 3 



1 



!(£) 

3(2 -(D) 3 
if £ > 2. 



+ !(£) 3 ifO<r<l 
if 1< £ < 2 



(17) 



As usual, in the Euler equation l|12|) an additional vis- 
cosity term is introduced to describe more complicated 



if VijTij > 



(19) 



with the following meaning of the symbols: Vij — Vj — 
VifCij — {cj + Cj)/2 is the mean sound velocity; hij = 
(hj+hi)/2,~p i j = (pj+pi)/2, a, (5 and r\ are constant coeffi- 
cients with typical values ~ 1,^2, and ~ 0.1, respectively 
(the latter coefficient is inserted only to prevent numerical 
divergences) ; and /y repre sents a shear corr ecting term to 
the artificial viscosity fcfr. ISteinmet dfl996|) . 

The code calculates the time-step of each individual 
gas particle by means of the Courant condition, 



At 



Chi 



C.i 



[hi\Wvi\ + Cj + 1.2(acj -I- /3 m axj\pij\)] 



(20) 



where C is a constant nearly equal to ~ 0.3. For all kinds 
of particle (both collisional and collisionless) in presence of 
gravi ty a more stringent condition is required ijKatz et alJ 
Il99ri) 



Atc.i — rj x min 



1/2- 



(21) 



where e is the softening parameter and rj is constant with 
value ~ 0.5. Finally, the time-step for the i-th particle is 
the least of the two: At{ = min(Atc,i, Atc.i)', and the time 
step for the whole system is taken to be the smallest of all 
particles time steps. 

Furthermore, the smoothing length is v aried with tim e 
and space according to the prescription bv lBend 1^.990) 



dh 1,_ 
— = --hVv. 
dt 3 



(22) 



The softening parameters we have adopted for the differ- 
ent components of our simulations are given in Table 

3.1. Cooling 

Gas cooling by many radiative processes perhaps plays the 
dominant role in the collapse of baryonic matter in galax- 
ies and star formation, since it is the comparison between 
the cooling and the dynamical time scales, respectively 
given by 



_ E 

^cool ~ r ~ 



ZpkT 



E 2pA(T) 



7T / 2GM 
and t dyn ~-^—^- 



(where p is the mean gas density and A(T) is the cool- 
ing rate of the gas at the temperature T) that decides 
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whether the gas will collapse nearly on the free-fall time 
(if tcooi < tdyn) or just contract through subsequent equi- 
librium conditions (if t coo \ > t^y n ). 

Many cooling processes are known. In general, the ra- 
diative cooling depends on density, chemical composition 
and temperature (from this latter very strongly). For T > 
10 i K the most efficient process is bremsstrahlung emission 
from the ionised plasma. In the range 100 < T < 1Q 4 K 
collision between H2 molecule and/or H atoms dominate 
via rotational and vibrational energy decay. Finally, for 
T < 10QK, t he dominant contri bution comes from CO 
molecules. See lChiosi et alJ l)l998|) for more details on the 
subject. 

Another important cooling mechanism is the inverse 
Compton effect, which is particularly important at high 
redshifts. The cooling rate of a gas with electron number 
density n e and temperature T, in a field of radiation with 
density p R and temperature Tr, is 



A 



Aa T n e p R (T - T R ) 



comp 



(23) 



Estimating the cooling time for matter immersed in the 
CMB we obtain 

_ [3w p m e (l + z)- 4 } 

z comp — fo „ x \^) 

{ap,0-T\lRPcrit) 

which is nearly equal to the dynamical time for col- 
lapse red-shifts greater than 7, independently of the 
galaxy mass. As our simulations begin at high redshift, 
the inverse Compton effect has been taken into ac- 
count. I n this paper we hav e adopted the cooling rates 
given bv IChiosi et alJ l)l998l and references therein) who 
have amalgamated results from different sources and 
to whom the reader should refer for al l details. The 

:n ad opted by Carraro et alJ 



(1998), 


Buonomo et al. 


Chiosi & Carrard 


(2002) 



tion and evolution. 



3.2. Star formation rate and initial mass function 

Star formation is among the most complicate and poorly 
known physical aspects of astrophysics despite its rele- 
vance in many issues. From simple minded arguments, 
there are at least three prerequisites for a gas (molecu- 
lar cloud) to be eligible to star formation: the gas has 
to be in convergent motion, i.e. the velocity divergence 
must be negative; the gas must be gravitationally unsta- 



ble, i.e. it must satisfy the Jeans condition r, 



sound 



> 



(where T soun d is the local sound velocity); the gas must 
be cooling, i.e. it has to verify the relation r coo ; -C Tff 
(see Section l3~TJ) . Then, the rat e of star f ormat ion (SFR) 
is customarily expressed by the ISchmidtl l|l959ft law 



dp* 
dt„ 



dfg 

dt n 



= c' P f (25) 

where c* is the so-called dimensionless efficiency of star 
formation, and t g is a characteristic time scale. Normally, 



SPH codes treat star formation simply implementing the 
Schmidt law in the computational language and trans- 
forming part of gaseous particles which satisfy the three 
conditions above in new, collisionless particles of differ- 
ent mass ("stars"), thus giving rise to a huge increase in 
the total number of particles. In order to avoid this cre- 
ation of new p articles w i th dif ferent mass as a result of 
star formation, iLia et all l|2002() proposed a new interpre- 
tation of the Schmidt law: the star formation rate given by 
eqn. I|25J1 was interpreted as the probability that at each 
time step a gas particle is instantaneously and fully turned 
into a star particle (thus losing its collisional properties). 
They showed that, averaged over a large fraction of the 
Hubble time, this probabilistic view of the star formation 
rate converges to the results of more conservative methods 
in which the masses of the newly formed stellar particles 
strictly follow the dictates of the chosen IMF law. The 
models we are going to present are based on this view of 
the star forming process. 

Furthermore, even if iBiionomo et alJ l)200C|) argued 
that the condition on the velocity divergence of gas parti- 
cles (negative) is not strictly required, the same arguments 
do not apply to our case because of the initial expansion 
phase of our proto-galaxies. Neglecting this, conditions, 
star formatio n could occur too early on. 

Following iBuonomo et all I 2000n . we adopt c*=1.0 for 
all the models. The characteristic time scale is chosen 
to be the maximum between r coo / and Tff time-scales. 
Numerical experiments show that in most situations t g = 
Tff is a good choice. 

IChiosi Carrarol §002) have found that the star for- 
mation history (rate vs time) depends on the depth of 
the gravitational potential well of a galaxy. In the case of 
deep potentials (such as in massive and/or dense galaxies), 
once star formation has started and energy is injected to 
the gas by supernova explosions, stellar wind, etc., this 
is not enough to push the gas out of the potential well. A 
sort of balance between cooling and heating is reached and 
the gas consumption by star formation goes to completion. 
Star formation cannot stop until the remaining gas is so 
little that any further energy injection will eventually heat 
it up to such high energies (temperatures) that the gravi- 
tational potential is overwhelmed. No more gas is left over 
and star formation is quenched. The star formation his- 
tory resembles a strong unique burst of activity, a sort of 
monolithic star forming event, taking place over a certain 
amount of time, of the order of 1 to 2 Gyr. In contrast, 
in a galaxy of low mass and/or density and hence shallow 
gravitational potential, even a small star forming activity 
will heat up the gas above the potential well. Some of it 
is soon lost in galactic wind, the remaining one becomes 
so hot that it will take long time to cool down and to 
form new stars. The cycle goes on many times in a sort 
of repeated bursting mode of star formation taking place 
during long periods of times if not for ever. 

In each generation of stars, otherwise known as Single 
Stellar Population (SSP), the stellar masses are known to 
obey a distribution function, named Initial Mass Function 
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Table 1. Initial parameters for Models A and B in the 
Standard CMD cosmological background 



Model 


A 


B 


Initial redshift 


50 


53 




1 


1 










H = 50 kmAlpc-'s' 1 


50 


50 


C"8 


0.5 


0.5 


Spin parameter 


0.02 


0.02 


Initial mean overdensity 


0.12 


0.12 


IMF 


Kroupa 


Kroupa 



Table 2. Initial dynamical and computational parameters for 
Models A and B. 



Model 


A 


B 


Initial number of gas particles 


13719 


13904 


Initial number of CDM particles 


13685 


13776 


Initial total mass (10 iz M©) 


1.62 


0.03 


Initial baryonic mass fraction 


0.10 


0.10 


Initial radius (kpc) 


33 


9 


Softening parameter for gas (kpc) 


1 


0.5 


Softening parameter for DM (kpc) 


2 


1 



(IMF), which is usually represented by a power law. The 
most popular of these functions is the one by ISaloeterl 
l)l955|) . More recent ones an d currently used in galaxy evo- 
lution models are those bv lArimoto fc Yoshiil l|l987|) and 
Kroupal l)l998|) . The three IMFs are shortly summarized 
here: 



$s(M) = CsM- 1 - 35 Salpeter 

( C K iM-°- 5 if M < 0.5 
<S> K {M) = \ CkM' 1 - 2 if0.5<M<l Kroupa 
[CkM- 1 - 7 ifM>l 

$a(M) = C A M- im Arimoto & Yoshii (26) 

where C<?=0.1716, Cxi=0.48, C K =0.295 and C A =0.145. 
Needless to say, the IMF bears very much on many aspects 
of galac t ic evo lution. In our models we have adopted the 
IKroupal l)l998(l IMF, originally designed to describe the 
situation for the solar v icinity. We prefer not to use the 
lArimoto & Yoshiil (|l987h IMF, even if it was claimed to 
be best suited to early- type galaxies, because it would 
produce a too high metallicity (see below). 

3.3. Heating 

Gas heating is due to supernova explosions, stellar winds 
and UV radiation emitted by massive stars, cosmic ultra- 
violet light, and to the cosmic background radiation. In 
addition to this we have several mechanisms of mechani- 
cal and dynamical nature which are already accounted for 
by the first term of the energy in eqn. . 



The total heating rate by radiative processes is given 

by 

r _ E S Ni + Esnii + E w + E uv 

At 1 ' 

where the terms refer to Type la and Type II super- 
novae, stellar winds and ultraviolet flux, respectively. In 
our model we consider only the supernova heating and 
neglect the other two. In brief, (i) the UV flux from mas- 
sive stars is absorbed by the interstellar molecular gas and 
re-emitted in the far infrared and thus is soon lost by radi- 
ation; (ii) the kinetic contribution by stellar wind should 
be treated on much smaller scales than those permitted by 
the resolution of these simulations. This is a point to keep 
in mind because the energy budget from stellar winds from 
massive stars may parallel that by supernova explosions. 
Finally, we do not include the heating by UV cosmic radia- 
tion. This may be a point of weakness because this heating 
source could delay gas cooling and star formation in turn, 
and even stop star formation at late r times by contrasting 
gas co oling. Although according to iNavarro fc Steinmet d 
it can be neglected at least during the first stages 
of galaxy formation, careful investigation of this issue is 
mandatory to properly understand the mechanisms of star 
formation. 

Supernova explosions return nuclearly processed gas 
to the interstellar medium, and heat it up. Type la super- 
novae from binary systems and Type II supernovae from 
massive stars affect in a different way the heating of the 
interstellar medium because of their much different evolu- 
tionary time scales. The progenitors of Type la are CO- 
white dwarf accreting mass from a companion, which in 
turn can be either a red giant or another CO-white dwarf. 
When and if the mass accreting white dwarf grows to the 
Chandrasekhar limit 1.4 Mq, C- ignition can take place 
inducing thermal runaway and consequent supernova ex- 
plosion. The time scale for this to occur ranges from a 
few hundreds Myr to several Gyr. In contrast, Type II su- 
pernovae generated by core collapse in massive stars have 
time scales of the order of a few Myr. Taking into ac- 
count such different behaviour is an important aspect of 
modelling galaxy evolution. We assume that each super- 
nova explosion liberates k$N x 10 51 erg of energy; anyway, 
only a small fraction of w hich - say 10% - can be thought 
as given back to the ISM ^Thornton et alJ ll998): the pa- 
rameter kg n fixes the amount of energy which is actually 
made available to the gas and not dispersed by radiative 
processes. This energy injection increases the internal en- 
ergy of the gas particles in the proximity of the explosion 
according to the SPH formalism and some kinematical / 
dynamical effects ought to be expected. However, as for 
stellar winds, it is difficult to correctly take into account 
the effect of this on the kinematics of the nearby gas parti- 
cles, as the scales on which these effects should be notice- 
able are too small if compared to the resolution of our sim- 
ulations. Therefore we simply give back the whole amount 
of energy to the thermal budget of the nearby particles. 
This approximation has little impact on the dynamical 



E. Merlin & C. Chiosi: Formation and evolution of early-type galaxies 



9 




Fig. 1. Positions on the xy plane of the Dark Matter par- 
ticles for Model A (proper Mpc) at different times; from 
left to right and top to bottom, z = 7.6, 3.2, 1.6, 0. 
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Fig. 3. The same as in Fig^but for the star particles. 



the chemical yields per stellar generation, second to suit- 
ably sprea d the newly form ed metals into the surround- 
ing gas; see lLia et al.l (j2002^ for a complete description of 
the physical processes involved in the computation and of 
their implementation in the code (the laws governing the 
gas restitution together with the chemical enrichment arc 
treated in the same probabilistic manner as the ones gov- 
erning the star formation and supernova rate processes) . 

The metals are distributed among the gas particles by 
means of the diffusion equation 



dZ 
~dt 



-k\7 2 Z 



(28) 



Fig. 2. The same as in Fig. Q]but for the gas particles. 



where n is the diffusion coefficient, for which we adopt 
the estimate bvlLia et alJ ll2002l) k = 9.25 x 10 16 km 2 s" 1 
based on the iThorntonet al.l l)l998|) models. In our sim- 
ulations we explicitly follow the contribution to the mass 
abundance often elements, namely Z (metals in general), 
He, C, O, N, Mg, Si, S, Ca and Fe. 



evolution of the whole system, due to the high cooling 
rates involved in the pro cess (see e.g. lCarraro et al.H l998: 
iKawatafc CihsonlfcOOlfl . 

To describe the rates of supernova explosions we a dopt 
the popular formulation by Greggio & Renzini ( 1983}) ; see 
there and references therein for details. To describe the 
energy feed back, we adopt on ce more the probabilistic de- 
scription of Lia et al to which we remand for a 
detailed description. 

3.4. Chemical enrichment 

The material, processed by nuclear reactions in stellar in- 
teriors and given back to the interstellar medium by su- 
pernova explosions and/or stellar winds, enriches the gas 
in metals. To correctly estimate the chemical enrichment, 
first one has to carefully calculate the stellar ejecta and 



4. Galaxy models: the parameters 

4.1. The cosmological initial conditions 

As already pointed out, in this paper we present two sim- 
ulations, whose parameters are summarized in Tabled 

All simulations are framed in a flat Universe with spec- 
tral index for the initial perturbations n=l. Denoting with 
fii the ratio between density of the generic component 
i and the critical density to close the Universe, in our 
Standard CDM models we adopt f2& = 0.1 (for baryons) 
and f2 c = 0.9 (for Dark Matter). The Hubble constant is 
taken to be 50 km sec -1 Mpc -1 . 

Little is said in literature about the choice of the ini- 
tial temperature to be assigned to baryonic particles. In 
general, the choice depends on the temperature coverage 
of the cooling processes under consideration, and very of- 
ten the initial temperature is set at the lowest value of the 
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Fig. 4. Position on the xy plane of the Dark Matter par- 
ticles for Model B (proper Mpc) at different times; from 
left to right and top to bottom, z = 6.3, 4.2, 2.2, 1.0. 
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Fig. 5. The same as in FigQJbut for the gas particles. 

temperature interval. This is less of a problem with our 
code becaus e we have implemen ted the large tabulations 
supplied by IChiosi et al.l (|l998fl which extend down to 
molecular processes at low temperatures (as low as about 
10 - 100 K). 

A rough estimate of the mean background temperature 
of the Universe at the redshift we are interested is derived 
as follows. At the decoupling redshift z ~ 1000, the tem- 
perature is about 3000 K. After decoupling, the Thomson 
scattering ceases, but a sufficiently large fraction of free 
electrons remains to keep matter in equilibr ium with rad i- 
ation via Compton scattering until z ~ 160 l)Madail2003h . 
Since then the temperature ever decreases (as far as mat- 
ter component are concerned) proportionally to (1 + z) 2 . 
Given that we are considering over-dense regions, we as- 
sume the initial temperature of the sphere to be equal 
to the mean temperature of the Universe, at the time in 



Fig. 6. The same as in Fig0]but for the star particles. 

which the mean background density was equal to the mean 
density of our proto-galactic sphere. With the initial red- 
shifts we have adopted (in the range 50 < z < 100), the 
temperature falls in the interval 10 to 150 K. 

4.2. Final parameters for the galaxy models 

For all simulations the initial smoothing length parameter 
for the baryonic component has been set equal to ten times 
the softening length. The softening length e in turn has 
bee n chosen in su ch away that both a good resolution and 
the lEvrardl l)l98S|) relation are secured. This latter requires 



Gi 



GM 



(29) 



i.e. e ^> R/Npart- In Tableware listed the dynamical and 
computational parameter we have adopted for the three 
models (the cosmological parameters have already been 
given in Table 

The basic units for length, time ad mass used to ex- 
press all physical quantities entering the models are: 

— length: 1 proper Mpc; 

— time: the lifetime of the Universe, corresponding to the 
cosmological model adopted for the simulation; 

— mass: the unitary mass is obtained supposing the free- 
fall time-scale of a sphere with unitary radius equal to 
the unit time so that the gravitational constant G = 1 
in code units. 



5. Galaxy models: the results 

Model A has been evolved up to age of about 13 Gyr (i.e. 
the age of the Universe), whereas Model B has been fol- 
lowed up to 5 Gyr. Nevertheless, the last computed model 
is well beyond the stage of maximum stellar activity, i.e. 
the galaxies have already settled down to a passive evolu- 
tionary regime. 
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Fig. 7. Top Panel: Final surface density profile for 
Model A (log[g/cm 2 ] vs kpc). The dots are averaged sur- 
face densities projected on the xy plane. The solid line is 
a Sersic law (see text for details). Bottom Panel: The 
same but for Model B. 

5.1. Spatial evolution 

The series of Figs. ^ through show the positions of Dark 
Matter, gas and star particles projected on the xy plane 
at different values of the rcdshift for Models A and B as 
indicated. It must be noticed that the panels for Model A 
has length scale of 300 kpc along both axes, whereas the 
length is 100 kpc in the case of Model B. 

One may sketch the following evolutionary picture for 
the two galaxies. The proto-galaxy initially expands with 
the Hubble flow and the expansion becomes slower and 
slower until it detaches from it first in the central regions 
(or in general in the denser regions) and later in the out- 
skirts. Not all gas is able to collapse but a certain (often 
huge) fraction remains at the edges of the just born galaxy. 
This is the result of the inhomogeneous distribution of the 
initial densities that lead some regions to collapse earlier 
than others independently of their distance from the cen- 
tre of the system. As a whole, the collapse deviates from 
the ideal spherical case. In particular, one may stress that 
the final baryonic-to-Dark ratio will be different from the 
initial one, in the virialized region of the galaxy. In Model 
A the density peak is nearly uniformly distributed around 




0.5 1 1.5 2 2.5 

log r 




0.5 1 1.5 

log r 



Fig. 8. Top Panel: Final spherical density profile for 
Model A (log[g/cm 3 ] vs log[kpc]). Dots and squares are 
respectively CDM and stars averaged densities in spher- 
ical shells at the corresponding radius. The solid line is 
the Hcrnquist law; the dotted line is the NFW profile. 
Bottom Panel: The same but for Model B. 

the geometrical centre; it forms a massive peak and sev- 
eral small sub-units. In Model B we note the evolution of 
two distinct units, one evidently smaller than the other, 
which eventually merge to form a single object; the merger 
is completed at z ~ 1.3. 

Once the collapse has started, while the collisionless 
Dark Matter begins to form clumps of high density the 
gas is decelerated and heated up by its own viscous fric- 
tion which causes the onset of shock processes; the kinetic 
energy of the gas is largely turned into internal energy by 
such shocks; strong cooling begins to occur, thus initiating 
star formation. Baryonic matter, in the form of stars, can 
now freely fall in the potential wells of Dark Matter. 

The first stars, in general in small isolated groups, start 
to form at z ~ 12 in Model A and z ~ 15 in Model B. 
Star formation proceeds and as shown below it occurs in 
a sort of broad burst whose duration is of the order of one 
or two Gyrs; then it slows down until it (almost) stops. 

Following the overall collapse, the energy injection by 
shocks and stellar feedback increases the energy (tem- 
perature) of the gas left over by stars to the point that 
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Fig. 9. Top Panel: Star formation rate (-M© per year) in 
Model A. Bottom Panel: The same but for Model B. 

part of it can overwhelm the gravitational potential. The 
gas escapes the galaxy and star formation is eventually 
quenched off. 

At the end of the galaxy and star formation period, 
the morphology of all our models seems to nicely resemble 
that of real galaxies. 

To highlight the structure of our model galaxies, first 
we analyse their three-dimensional shape. To this aim we 
chose a suitable cut-off radius and consider the stellar and 
Dark components inside this radius. The axial ratios can 
be roughly estimated as 



! Zmizf 



(30) 



where x, y, z are the coordinates along the axis of the re- 
sulting ellipsoid. For Model A, inside the half mass ra- 
dius, we find the axial ratios b/a — 1.08 (stars) and 1.14 
(Dark Matter), and c/a = 1.07 (stars) and 1.17 (Dark 
Matter). For Model B we get b/a = 1.04 (stars) and 1.14 
(Dark Matter), and c/a = 1.00 (stars) and 0.96 (Dark 
Matter). We can therefore state that, in general, the stars 
are nearly spherically distributed, while the Dark Matter 
haloes show a weak asymmetry. Owing to the smallness of 
the deviation, the two models arc considered as spherically 
symmetric. 



Fig. 10. Top Panel: Formation time (log[Gyr]) vs radial 
distance (log[Mpc]) of the stars in Model A. It is easy to 
notice how the galaxy we see today is gradually built up 
moving from the external regions toward the centre where 
star formation lasts very long. Bottom Panel: the same 
but for Model B. 

The panels of Fig. show the final surface density 
profiles of the stellar contents of Models A and B projected 
onto th e xy p l ane w hen observed face on. The solid line 
are the ISersid l|l968j) surface density profiles 



£ = S e exp{- (0.324 



2m)[(R/R eB ) l ' r > 



1]} (31) 



where R e B, the effective radius, is the radius containing 
half mass of the galaxy (in this case referred to cylindri- 
cal symmetry), S e is the mass density inside this radius, 
and to is a numerical parameter which correlates with the 
galaxy absolute magnitude; for m = 1 one finds the ex- 
ponentia l law of dwarf elliptic als, while for to = 4 one 
gets the Ide Vaucouleurd l)l948|) law. The agreement be- 
tween the numerical results and the de Vaucouleurs law is 
remarkably good for Model A; only in the very central re- 
gions (i.e. r < 2 kpc) the density of star particles seems to 
be too peaked. Model B seems well fitted by a Sersic law 
with to ~ 6. The best- fits have been made minimizing the 
dispersion of the points \/(pi{r) — Pfit(r)) 2 ) with respect 
to the analytical relationship. As noticed by ICaon et alJ 
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Fig. 11. Variation of (from top to bottom) kinetic, ther- 
mal, total and potential energy in code units as function 
of time in Gyrs for Model A (all components) . 





Fig. 13. Top Panel Type la supernova rate (explo- 
sions per year) as function of time for Model A. Bottom 
panel: the same as above but for Type II supernovae. 



Fig. 12. Top Panel: Trend of the Virial trace 
— 2Ekin/ E pot for Models A and B from top to bottom. 
Ages are in Gyr. 



(1993), one would actually expect smaller values of m for 
smaller galaxies, while we obtain the opposite. This is a 
drawback of the present models, which however should not 
be rejected because of it. Profiles with m=6 and m=4 are 
indeed quite similar, so that most likely the present sim- 



ulations do not possess the resolving power to distinguish 
between the two. It is safe to say that the general trends 
of the mass profiles resemble those of real galaxies, but 
a detailed comparison between theory and observations is 
not yet possible with the present models. 

Considering the models as roughly spherical, we try to 
fit their three-dimensi onal density profil es with theoretical 
analytical laws, i.e. the lHernauistl 1 1990) law for the stellar 
component 



p(r) 



Ma 



1 



(32) 



2tt r (r + a) 3 

where M is the total mass of the stellar system and a 
is a scale length that de pends on the half-m ass radius 
r i/2 = (1 + \/2)oj and the lNavarro et alJ l)l996j) universal 
profile for the Dark Matter component 



p(r) 



(33) 



Pent (r/r s )(l + r/r s ) 2 

where p cr it is the critical density, r s is a scale radius at 
which the slope of the density profile 133|) is —2, and S c 
is a characteristic (dimensionless) de nsity that depend on 
the cosmological model and r s (see iNavarro et all Il99fil 
for details). 
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Fig. 14. Top Panel: Chemical evolution of the star and 
gas particles in Model A as a function of the age; left 
panels show the mean total metallicity Z in solar units, 
right panels the ratio [O/Fe]. The age is in Gyrs. Bottom 
Panel: The same as above but for Model B. 

Fig.Elshows the mass densities for our stellar and Dark 
Matter components compared with the theoretical profiles 
(fits are obtained the same way as above): the agreement 
seems quite good in the middle regions. Obviously, little 
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Fig. 15. Number distribution of stars in metallicity at the 
end of the simulation for Model A. 
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Fig. 16. Top Panel: Metallicity (log[Z/Z Q ]) vs age of 
formation (Gyrs) of the stars in Model A. Bottom Panel: 
The same but for Model B. 

can be said about the most central regions (i.e. less than 
~ 2.5 kpc), due to the effects of the so ftening length and, 
in general, of numerical nature (see e.g. lPower et al.lf2003l 
for a discussion of these issues). 

For both models the Dark Matter profile is very 
steep in the innermost central regions. Current cosmolog- 
ical si mulations give profile s with inner sl opes between 
~ - 1 ((Navarro et alJ Il99fij) and ~ -1.4 l|Moore et alJ 
199§j). However, there are evidences that the inclusion 
of non-adiabatic treatment of the gaseous component, 
and particulary the star formation pro cess, could steepen 
the Dark Matter density profile (e.g. | Lewis et all 12 000: 
iKawata fc GibsoiJl2003l bnedin et "al1l2004^ . IPower et alJ 
1)200.^ also argue that poor temporal resolutions could 
give rise to artificial cusps in the innermost regions. 
Therefore, considering the uncertainty still affecting the 
whole issue, the present results are acceptable. 

Furthermore, we note that, contrary to common expec- 
tations, in the central regions of our models the density 
of Dark Matter exceeds that of baryonic matter (stars): 
the Dark Matter halo seems to be strongly concentrated. 
We ad opt the definition of concentration bv lNavarro et alJ 
(1996) as the ratio c = R V i r /r s between the virial radius 
R V i r of the Dark Matter halo and the scale radius of rela- 
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Fig. 17. Radial gradient in the spherically averaged 
metallicity for the stars of Model A (top) and B (bottom) . 
Radii are in log[Mpc]). 



Fig. 18. Top Panel: Gas density (log[g/cm 3 ]) vs tem- 
perature (log[K]) of Model A. Mid Panel: The same but 
for Model B. 



tion i)-'S-SI) r s . As virial radius we adopt the radial distance 
at which the mean total density of the systems is 200 
times the background density; so, in this estimate both 
Dark and Baryonic matter are included. However, owing 
to the small fraction of baryons with respect to the total 
mass, the results does not significantly depend on their 
presence. The scale radius is derived from the best-fit of 
our numerical results to the NFW profile (see Fig. . The 
concentrations we find for Models A and B are ~ 50 and 
~ 20, respectively. 

These values are higher than current expectations 
which indicate conce ntration of ~ 10 for galactic haloes 
I Navarro et alJll99ft >. This drawback of our models re- 
quires a careful future analysis. First of all, as already 
mentioned, the presence of stars can steepen the Dark 
Matter profiles. In this contest the effects of different 
star formation recipes and/or different choices for the 
model parameters should be investigated. Second, quasi- 
cosmological initial conditions ignore the possible late- 
time infall of gaseous matter from the IGM into the virial- 
ized halo of the new born galaxy. Moreover, the observa- 
tional estimate of the mass profiles in the centre of galaxies 
is derived with the aid of (theoretical) mass-to-light ra- 
tio for the stellar components, which depend on the IMF 



as well as stellar evolution properties. While these lat- 
ter are somewhat established, this is not the case of the 
IMF. Passing from one IMF to another m ay easily change 
the m ass-to light ratios by a large factor l|Portinari et alJ 
2004J). In addition to this, our galaxies seem to be com- 
pletely relaxed and in equilibrium (see below). This prob- 
ably means that a baryon-dominated centre is not strictly 
necessary to secure the stability of the systems. Last, the 
density profiles of the lumi nous component of our mod- 
els fairly agree with those of lHernauistl l)l990|) . Therefore, 
as far as we can tell, neither theoretical nor observational 
evidences actually exist constraining the central regions 
of galaxies to be always baryon-dominated, each galaxy 
being peculiar in its formation and virialization history. 
In conclusion, we feel confident to consider our models as 
reliable representations of the gross scale features of real 
galaxies and the grand design of galaxy formation in which 
these numerical simulations are framed. 

Finally, the effective radius of the galaxy is calculated 
using the Hernquist density profile, so that 

^ eB — ~ToT' ( 34 ) 
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Fig. 19. Top panel: Radial velocity (km/sec, points) and 
escape velocity (dotted line), vs radius (Mpc) for Model 
A. Both velocities are in proper km/s. Bottom panel: 
The same but for model B. 

Table 3. Masses and half mass radii for the various 
components of Models A and B. Masses are in units of 
10 12 M Q , radii are in kpc. 



Model 


A 


B 


Final star mass M s 


0.091 


0.0029 


Final gas mass (collapsed) 


0.062 


0.0004 


M a /M bar 


0.56 


0.82 


Half mass radius of stars 


7 


1 


Effective radius of stars 


5.2 


0.8 


Sersic law index m 


4 


6 


Half mass radius of DM 


52 


15 


Concentration (NFW profile) 


50 


20 


Virial radius 


300 


41 


Age of the last model (Gyr) 


13 


5 



toernauistlll99fl . Table H gives the values of the parame- 
ters concerning the global properties of the density profiles 
of the models. 

5.2. The Star Formation History 

The history of the star formation rate, in solar masses 
per year, is shown in Fig. El for Model A (top panel) 



and B (bottom pa nel). They can be compar ed, f or exam- 
ple, with those by IChiosi fc Carraro! {2002) and iKawatal 
l)l999l EoOlalbh . It is evident that in all cases the galaxy 
goes through a phase of strong stellar activity, in a sort 
of burst-like mode lasting about 2 Gyr, followed by very 
slow activity (of the order of 1-5 M /yr) that continues 
for long periods of time. Remarkably the bulk of stars are 
in place at the galaxy age of 2.5 Gyr, i.e. at the redshift 
z = 1.95 for Model A and z = 2.2 for Model B. 

This behaviour shares the properties of the hierarchical 
and monolithic scheme. It is hierarchical because in the 
proto-galaxy lumps of matter merge together to assemble 
the galaxy we see today. The merging history is complete 
at redshift z = 2.5 for the Model A and z = 2.2 for Model 
B, except one major merger which happens at z ~ 1.3. It 
is monolithic because the largest fraction of the star mass 
is built up very early before z=2, for all models. 

The results we have obtained clearly favour the mono- 
lithic scheme of galaxy formation (cfr. Section [T]) as a 
unique episode of star formation during very early stages 
can generate a structure very similar to what we see to- 
day; even in the case of a subsequent major merger as that 
of Model B, the star formation history does not undergo 
substantial changes, perhaps due to the small amount of 
gas left over by the first burst of activity. 

How star formation actually works obviously depends 
on the initial conditions. In our models, the results favour 
the mergers of small sub-units in very early epochs rather 
than the formation of big galaxies from mergers of other 
galaxies at much later times. 

A huge fraction of the collapsed gas is turned into 
stars and also a considerable amount of gaseous material is 
blown away before it can collapse in the proto-galaxy po- 
tential well. With the present recipes, the efficiency of gas 
restitution by dying stars refuels the interstellar medium 
with sufficient new gas to engender minimal star formation 
after the galaxy has reached fully relaxed conditions. This 
is less of a problem here because the poor mass-resolution 
makes the effects on the star formation rate of the trans- 
formation of one or two particles per step more dramatic 
than they actually are. Moreover, slightly different condi- 
tions on the recipe for the star formation rate could get rid 
of it. In relation to this, many experiments have recently 
been made to evaluate and simulate the effects of stellar 
energe tic feedback on the interstel l ar medium (cfr. for in - 
stance ISpringel fc Hernauistll2003t iMarri fc White! l2003h . 
Sufficiently strong feedback would stop the star forma- 
tion activity by evaporating the cold gas clouds in which 
new stars would form. However the poor resolution of the 
present numerical simulations does not allow us to prop- 
erly investigate this problem. As already pointed out, we 
are not able to reproduce the effects of stellar re-heating 
on the kinetics of gas particles, so that we must limit 
ourselves to just fuel the thermal budget, which actually 
proves to be ineffective, due to the strong cooling phe- 
nomena which dominate the gaseous non-adiabatic pro- 
cesses at the temperatures we are dealing with. Therefore, 
past the strong, dominant episode of star formation, the 
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long tail of stellar activity can be simply ignored for all 
practical purposes, for instance the derivation of spectro- 
photometric properties. 

As already mentioned, the by far dominant episode of 
star formation occurred at very early epochs and lasted 
for a significant fraction of the Hubble time. The duration 
of it (about 2 Gyr) is longer in Model B than in Model A 
thus leading to a higher ratio between the total baryonic 
mass and the left-over gas mass (cfr. Table [2J ■ AH this 
bears very much on the chemical evolution of the (model) 
galaxies as discussed in Section ET41 

Table |3 lists the mass of the various components at the 
end of the simulations together with the half mass radius 
of the stellar and Dark Matter distributions. CDM turns 
out to be strongly concentrated, as already noticed. The 
final stellar-to-gas ratios (considering only the collapsed 
and/or reprocessed gas) turns out to be ~ 1.47 and ~ 7.25 
for Models A and B, respectively; as already pointed out, 
Model B results much more efficient in forming stars, leav- 
ing only ~ 13% of the baryonic collapsed particles in the 
gaseous form (see Table [21 for details). The mass of Model 
B is not small enough to consider it a dwarf elliptical, but 
this results is still a bit surprising, since we would expect 
that a less deep potential well would produce less effi- 
cient star formation processes (see e.g. IChiosi fc Carrarol 
120021) . Remarkably, the vast majority of stars in Model B 
are formed before the major merger at z ~ 1.3, thus rul- 
ing out the possibility that the merger could be respon- 
sible for the different efficiency in forming stars (indeed, 
a small increase in star formation can be noticed in coin- 
cidence with the merger, but after the formation of the 
bulk stellar content stars). A possible explanation can be 
sought in the presence of two distinct potential wells. The 
central regions of an even small potential well can indeed 
be more efficient in forming stars, than the peripheral re- 
gions of a bigger one. In other words, while the gaseous 
particles in the outskirts of a single, deep potential well 
begin to heat up by stellar winds and supernova explo- 
sions before they can fall inside the well and therefore are 
likely subtracted from the star forming activity, in pres- 
ence of two (or more) potential wells a larger percentage 
of gas particles, in different regions of the forming galaxy, 
can be turned into stellar particles. 

FigllUI shows the age of the stars as a function of their 
radial position for Models A and B. It clearly illustrates 
how a galaxy is assembled. The star formation process in 
the central regions goes on for very long periods of time; 
younger and more metal-rich stars are therefore expected 
to be found in the core of galaxies, whereas older and 
more metal-poor stars are likely to be found in the outer 
regions. 

5.3. Energy 

The panels of Fig^] from top to bottom, show the vari- 
ation as a function of time of the kinetic, thermal, to- 
tal and potential energies (all expressed in the physical 



units we have adopted in the code) of Model A. The same 
quantities for Model B have similar behaviour. The trends 
resemble nicely what we expect to see in an expansion 
and subsequent collapse of an isolated system; notice that 
the thermal energy, which is provided only by the gas, is 
very small when compared to the kinetic and the potential 
terms. The total energy is not strictly conserved (only a 
few per cent), as the inclusion of non-adiabatic processes 
such as cooling causes the loss of part of the energetic 
budget of the gaseous component. 

Fig. ITS! shows the Virial Trace V — —2E Kin /E pot as 
a function of time for the two models. It is easy to follow 
the relaxation of the system, which tends to a state of 
equilibrium after some small oscillations (if the systems 
had been strictly collisionless and more symmetrical these 
oscillation would have been much stronger, as we expect 
in a violent relaxation process). The expected value for an 
ideal system would be V = 1, but here the presence of 
viscosity and cooling of the gas tend to shift the ratio to 
somewhat lower values to the loss of energy. It is worth 
noticing in Model B the bump in coincidence with the 
major merger occurring at ~ 3.5 Gyr. 

Moving now to heating processes, the rate of supernova 
explosions of Model A is shown in Figs. ED (the trend is 
nearly the same for Model B). The top and bottom panels 
are for Type IA and Type II supernovae, respectively. As 
expected the rate of Type II supernovae strictly follows 
the star formation rate: as massive stars are short lived, 
as soon as star formation starts and/or ceases the same 
does the Type II supernova rate. The case of Type la is 
different because they require a certain amount of time to 
show up and continue to explode for very long time after 
the end of the star forming period. 



5.4. Metallicity 

Fig. El displays the evolution of the mean metallicity Z 
and mass abundance of Fe and O both for stars and gas 
in Models A and B; the ratio [O/Fe] is in the standard 
spectroscopic no tation. The plots c an be compared with 
those in Fig.5 of lChiosi fc Carrarol l)2002|) . 

Looking at the [O/Fe] vs age relation displayed in Fig. 
1141 we notice that the so-called a- enhancement problem is 
present. Both models have the stellar content strongly en- 
hanced in a-elements, however the degree of enhancement 
is lower in Model B (lower total mass) and in both cases 

decreases with the age. These trends are ag ree with the ob - 

I — . — I J Ti 

servational information, see the review by IC hiosi (2000), 
and the s imilar analysis made bvfTantalo fc Chiosil 1)2002 ) 
using the lChiosi fc Carrarol l|2f)0^ models. 

The final mean stellar metallicity in Model A is nearly 
solar (< Z >~ 0.8Z©). Remarkably also Model B shows 
similar values of mean metallicity even if its evolution was 
stopped much before the present age, thus showing how 
the chemical enrichment seems to be almost completed in 
the first Gyr of a galaxy life. 
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Fig. 1 151 shows the final number distribution of stars per 
metallicity bin (Z/Z Q ), limited to Model A: 4% of stars 
have mean total metallicity lower than 1/3 solar, 67% have 
mean total metallicity in the range 1/3 solar to solar, a 
significant fraction goes up to three times solar, and a tiny 
fraction even up to 5 times solar. The low fraction with 
metallicity lower than 1 /3 solar secures that the so-called 
G-dwarf problem does not occur, whereas the tiny frac- 
tion around 5 times s olar could exp l ain th e UV excess in 
elliptical galaxies. See lBressan et al 1 lll994l) for a thorough 
discussion of these topics. 

At any age a large scatter in the metallicities of the 
star particles has to be expected. This is shown in Fig. 1161 
The situation resemb les the one observed amon g the stars 
of the Galactic Disk ijEdvardsson et al.lll993bllali . 

Finally, in Fig. 1171 we display the spherically averaged 
radial gradient in metallicity for Model A and B; stars 
closer to the centre are also m ore metal rich. For elliptical 
galaxies, iDavies et al 1 lll993h find the metallicity gradient 
Alog{Z) I Alog{r) ~ -0.2 ± 0.1. The gradient of Model 
A is ~ —0.30, in good agreement with the observational 
estimate. Model B shows instead a steeper gradient, ~ 
—0.50, due to the larger amount of stars produced towards 
the centre. 

As a last remark we like to report here that similar sim- 
ulat ions with almost identica l initial conditions but with 
the lArimoto fc Yoshiil l)l987|) IMF yield very high values 
of the final mean metallicities, i.e. up to 5 or 6 times the 
solar abundance llMerlinl l2005h . The reason is that the 
lArimoto fe Yoshiil l|l987|) IMF favour massive stars and 
Type II supernovae and ther efore predict m uc h higher 
rates o f gas restitutio n than thelKroupal l|l99Sj) or lSalpeterl 
(|l955h IMF (see e.g. Ilia et al.ll2002l) . As such high mean 
metallicities do not agree with current observational es- 
timates, in the present study we prefer to abandon the 
lArimoto fc Yoshiil l|l987f) IMF even if it was originally tai- 
lored to model ellip tical galaxies. The rec ent calculation 
of stellar yields bv IPorti nari q\ ahl l)l998h which arc in- 
corporated in the lLia et al.l l|2002fl algorithm for chemical 
enrichment in N-Body PD-TSPH simulations, show that 
IMFs too skewed toward m assive s t ars ar e not required. 
The reader should refer to (Merlin! ( 120051) f or all details 
about these models with lArimoto fc Yoshiil {l987) IMF. 



5.5. Galactic Winds 

As already reported in Section ^ elliptical galaxies ought 
to suffer galactic winds in order to explain the colour- 
magnitude relation. To investigate this important issue, 
in Fig. 1181 we show the temperature vs density relation- 
ships for all gas particles of our models. The gas left 
over by star formation is made of four components - 
hot and cold, chemical ly enriched and unprocessed (cfr. 
IChiosi fc Carraroll2002|) . The hot gas, which is the domi- 
nant component, has relatively high densities: heated up 
by energy injection is about leaving the potential well of 
the galaxy. The fact the hot gas constitutes almost the 



totality of the gas content is due to the efficiency of the 
energy feed-back and the generally large total mass of the 
systems (relatively deep potential wells) which inhibits the 
occur r ence of earlier galac tic winds iGhiosi fc Garrard (cfr. 
l2002jk iKawata fc Gibson! (cfr. l2003h . Note that in Model 
B, as fewer gas particles are left over by star formation, 
the situation is less defined then in Model A, where we can 
clearly see the final dependence, for the hot component, 
of temperature on density. 

Fig. ^5] shows the radial velocity of the gas particles 
together with the radial profile of the escape velocity. The 
amount of gas (excluded the fraction which did not take 
part to the collapse of the proto-galaxy) which, at the end 
of our simulations, has the radial component of its velocity 
greater then the escape velocity at its distance from the 
centre of mass of the system, can freely leave the poten- 
tial well of the galaxy and can be considered as gaseous 
mass ejected as galactic wind. In all simulations there is 
a layer beyond which all the gas particles meet these re- 
quirement; the percentage of the gas which can be con- 
sidered as galactic wind is ~ 21% (Model A) and ~ 58% 
(Model B) on the total of collapsed gas. Note that these 
percentages are very different from one case to another; 
this could be due to several factors, among which the dif- 
ferent star formation history (as already pointed out, star 
formation is much more efficient in Model B than Model 
A, thus leaving a smaller amount of gas in the central 
regions of the system), and the different masses of the 
systems, since galactic winds are expected to be more effi- 
cient in galaxies of low mass llGhiosi fc Garra,roll200l . At 
large distances from the galaxies, some particles with al- 
most zero radial velocity can be found (see Fig. I19|) . They 
are the un-collapsed gas particles that did not take part to 
the galaxy formation process. They cannot be considered 
as galactic wind as they have not been blown away from 
the galaxies. 

Finally, we would like to note that a better treatment 
of the gaseous componen t should be con sider ed su ch as 
the multi-phase scheme of lMarri fc White] fe.g. l2003j) that 
would allow us to describe in a more realistic fashion both 
chemical enrichment and galactic winds. 



6. Conclusions 

In this study, with the aid of N-body simulations based on 
quasi-cosmological initial conditions, we have followed the 
formation and evolution of two early-type galaxies from 
the stage when they separate from global expansion of 
the Universe to their collapse to virialized structures, the 
formation of stars and subsequent nearly passive evolution 
up to the age of 13 Gyr. The cosmological background is 
the Standard CDM. We have highlighted the structural, 
dynamical and chemical properties of the models to be 
compared with observational data. 

The models contain several parameters, among which 
we recall 
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— The softening parameter for Dark Matter, gas and 
stars; 

— The specific efficiency driving the star formation 
rate; 

— The constant ksN defining the fraction of energy re- 
leased by a supernova explosion which goes into the 
gas as thermal energy. 

The uncertainty on these parameters opens the gate 
to future studies; better fine tuning of these quantities is 
required to get results in very close agreement with the 
observational data. Nevertheless the results we have ob- 
tained are satisfactory and in agreement with other studies 
of the same subject. 

Two systems made of Dark and baryonic matter, in 
current cosmological proportions, with positions and ve- 
locities for both type of particles derived from the spec- 
trum of cosmological perturbations suited to the cosmo- 
logical model under consideration have been investigated. 

All models started from the conform expansion of the 
Hubble flow, with the l owest rigid body rotation allowed 
by the CDM Universe ((Barnes fc EfstathiovJll987l) . They 
all have been followed for a large fraction of the Hubble 
time, through the phases of turn around and subsequent 
collapse of the galactic proto-clouds of gas to small clumps 
caused by the initial density perturbations to the for- 
mation of equilibrium structure made of Dark Matter, 
stars and remaining gas. Hydrodynamical and thermal 
processes caused by cooling, star formation, chemical en- 
richment, and heating have been included and followed in 
detail. 

The final systems morphologically resemble real ellip- 
tical galaxies of intermediate dimensions; this is proved 
by their su rface density pr o files m atching the theoreti- 
cal ones by Ide Vaucouleursl Jl948fl and/or ISersid lll968h 
that ar e expected for galaxies of the same mass and scale 
lengths. iHiorth fc Madsenl l(l99lh have indeed showed that 
violent relaxation in a deep pote ntial well leads to sys tems 
obeying the de Vaucouleurs law. lKatz fe Gunnl l(l99lf) and 
Kawat al lll999l) outline that dissipative processes, followed 
by the PD-TSPH code in great detail, bear very much on 
the final results. 

Star formation lasts longer in the central regions than 
in the outskirts thus producing gradients in metallicity 
that clos ely agree with the o bservational data for elliptical 
galaxies l|Davies et al.lll993|) . 

It is woth stressing that in any case the models still 
need to be improved in several aspects: (i) First the star 
formation that even if shaped in a huge initial burst fol- 
lowed by quasi-quiescence, never completely stops, be- 
cause of poor mass resolution and/or still inadequate 
prescription for the physical conditions at which stars 
can be formed. The long tail of minimal stellar activ- 
ity does indeed affect the spectro-photometric properties 
(e.g. colours) of the model galaxies that would disagree 
with observational data, (ii) Second, the density profiles 
that are not fully satisfactory, as in the central regions the 
baryons do not dominate the spatial density of matter as 



suggested by the observational data, (iii) Third, there's 
no right correspondence between the values of Sersic pa- 
rameter m to galaxies of different mass and the those de- 
rived from observational data, (iv) Finally, the predicted 
mass-metallicity relation as compared to the one inferred 
from observational data (colours, e.g. the CMR). The 
present model B with lower mass actually suffers more 
star formation and metal enrichment in turn. The rea- 
son for that is not yet clear. The present results cannot 
yet recover the m e tallici ty-mass relationship obtained by 
IChiosi fc Garrard ll2002h that fairly agreed with the ob- 
servational information. 

Nevertheless, it is worth noticing th at our main re - 
sults are in good agreement with those bv lKawatal 1 19991. 
who used the same method to generate the initial condi- 
tion s, as well as with those bv iKawata fc Gibson! l(2003h 
and iKobavashil ((2005) , who instead used the more refined 
method otherwise known as resimulation technique (see 
Section. l2.2() . This point is particularly relevant as it shows 
that our appr oach, developed i ndependently and in par- 
allel to that of IKobavashil (jioOfih . can be considered suffi- 
ciently accurate. 

Our simulations support the revised monolithic sce- 
nario of galaxy formation. Through the hierarchical col- 
lapse of Dark Matter haloes, clouds of gas meet, heat 
up, cool and eventually collapse forming stars in a single 
episode of star formation occurred in the far past, during 
the virilization of the proto-galaxy, perhaps followed by 
minimal stellar activity in the central regions and almost 
negligible morphological and structural evolution. 
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